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PREFACE 

This Memorandum is part of a continuing theoretical study of the 
geomagnetic field. The results should aid the development of repre- 
sentations of the geomagnetic field in space, and are directly appli- 
cable to both geophysical problems involving scalar and vector poten- 
tial fields, and to the design of conjugate point experiments. The 
work was supported by the National Aeronautics and Space Administration 
under Contract NASr-21(05) . 
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ABS TRACT 

The geomagnetic field is analyzed by spherical harmonics and by 
integrals. Series representations in spherical harmonics of geomagnetic 
field charts are compared for truncation of the series at the 6, 8, 

10, and 12 terms of degree. Scalings are at a uniformly spaced 
latitude-longitude grid from both U.S. and U.S.S.R. isomagnetic charts 
for 1955.0. A numerical integration method for analyzing the field 
is developed from Poisson's integral. A new surface grid, suitable 
for use with integral analysis, is described. This grid is based on 
subdivisions of a spherical icosahedron, and its points are almost 
uniformly spaced over a sphere. This integration method is applied to 
calculations of field values, field lines, and conjugate points. The 
results are compared with those of earlier spherical-harmonic analyses 
by Vestine and Sibley. A comparison is also made between those conjugate 
points calculated by spherical harmonics from different sets of 
coefficients derived from various sets of isomagnetic charts. Some 
minor but undesirable effects are mentioned that arise because of the 
uniform angular spacing of data points scaled from charts. The vari- 
ation of the earth's magnetic moment and the location of the dipole 
axis since 1835 is described and discussed. Finally, an extrapolation 

field into the earth's interior is described. 


of the geomagnetic 



-vii- 


ACKNOWLEDGMENTS 

The authors would like to acknowledge the assistance of 
Louise Kern and Anne Kahle in preparing the data and figures. 




-ix- 


CONTENTS 

PREFACE iii 

ABSTRACT v 

ACKNOWLEDGMENTS vii 

LIST OF FIGURES xi 

Section 

I. INTRODUCTION 1 

II. SPHERICAL HARMONICS 2 

III. INTEGRALS 4 

IV. SPHERICAL-HARMONIC ANALYSES OF ISOMAGNETIC CHARTS 

OF THE U.S. AND U.S.S.R. FOR 1955.0 9 

V. ANALYSIS BY NUMERICAL INTEGRATION 12 

Two Equal-Area Grids 12 

Numerical Formulae for Surface Integration 18 

Analytic Procedures 20 

VI. APPLICATION TO GEOMAGNETIC FIELD-LINE CAI£ULATION . . 26 

VII. COMPARISON OF CONJUGATE POINTS CALCULATED FRCM 

VARIOUS SPHERICAL-HARMONIC ANALYSES 31 

VIII . MAGNETIC MOMENT AND DIPOLE AXIS SINCE 1835 33 

IX. MAGNETIC FIELD IN THE EARTH'S MANTLE 35 

X. CONCLUSIONS 37 

APPENDIX: TABLES 39 

BIBLIOGRAPHY 63 



-xi- 


LIST OF FIGURES 


1 . 


2 . 


3. 


4. 


5. 


6 . 


7. 


8 . 


Method of generating an equal-area grid with latitude 

and longitude lines 

Method of generating an equal-area, equal-spacing grid 
based on the spherical icosahedron 


Coordinate system for surface integration over a 

single area sector (circular cap on sphere) 

Results of surface integral for circular cap at different 

positions relative to the cap 

Field lines calculated by means of spherical-harmonic 
and surface-integral analyses from 40°N , 90°E, altitude = 
0.1 earth radius. Step size along field line is 0-1 

earth radius 

Field lines calculated by means of spherical-harmonic 
and surface-integral analyses from 50°N , 0°E , 
altitude = 600 km. Step size along field line is 

0.25 earth radius 

(A) Geomagnetic north pole, from various analyses. 


1829-1960. 

25 -x/ 

(B) Moment, M, in units of 10 cgs 

Nondipole main-field components for 1955.0 at a depth 

of 1500 km. Contours are for downward (Z) component 

in gauss; horizontal component (H) is given by arrows . . 36 




- 1 - 


I . INTRODUCTION 

About 1839, Gauss first used spherical harmonics to analyze the 
geomagnetic field's potential function. Since then, at intervals of 
ten years or more, others have made such analyses (Chapman and Bartels . 
1940) . The results of their work, usually in the form of charted 
values, have seldom been rigorously comparable. This is hardly 
surprising, since they used different accuracies and distributions 
for the observational data points, assumed different numbers of 
spherical harmonic terms to fit the data, weighted their observations 
differently , and used different methods of analysis. 

This paper estimates the differences caused by defective data 
and procedures, and indicates the effect on a few major main-field 
parameters and their interpretation. It compares the computed results 
of spherical-harmonic analysis with those of integral analysis, in 
two meridional planes and at various heights above the earth's surface. 
It describes formulas and grids that are useful for integral analysis. 
It then compares the sample results with previous tabulations of the 
main field that have been used in analyzing particle data for the 
Van Allen radiation belts. It indicates the change since 1835 in the 
geomagnetic-pole position and in the earth's magnetic moment. Finally, 
it discusses extrapolation of the surface field into the earth's 


interior. 
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II . SPHERICAL HARMONICS 


Spherical harmonic analyses of the geomagnetic field usually 

represent the geomagnetic surface potential over a spherical earth 

(whose radius equals a) in this form: 

, . rH-1 


V = a £ — J (g™ cos mX + h™ sin mX) P™ (cos 0) (1) 

The earth’s center is taken as the origin of the three spherical 
coordinates: r, the distance from the origin; 9, the colatitude; and 
X, the longitude east of Greenwich (Chapman and Bartels , 1940). 
p™ (cos 0 ) are Schmidt’s semi-normalized associated Legendre polynomials 
of integral order m and degree n, and g™ and h™ are the Gaussian 
(Schmidt) coefficients. The north, east, and vertical (or downward) 
components of the surface magnetic field are then given by 


1 5V 15V 5V 

X = > X = — , Z = — . (2) 

r 50 • r sin 0 5X 5r 

The values of and h m are usually determined from the observed field 
to n n 

values. Most analyses provide values up to m = n = 6 and fit, by the 

o o 

method of least squares, weighted data taken from charts at 5 to 10 
intervals of latitude and longitude. An example is the analysis for 
1922 (Dyson and Furner , 1923); essentially the same methodology was 
followed by Vestine, et al . t for 1943, and for secular change at 
10-year intervals from 1912.5 to 1942.5 (Vestine , et al . , 1947). The 
results obtained are therefore comparable, since they include the 
influence of similar defects as well as advantages in methodology. 
Repeating the study of 1945 — but now using the data of Vestine, et al . , 
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up to m = n = 15 — gave surprisingly similar results (Fanselau 
and Kautzleben , 1958) . 

Other methods based on observed points (unequally distributed 
measurements at observatories) have afforded almost as good an 
approximation, though they have not well represented the field's 
distribution over the oceans. Analyses to terms of high degree 
(512 coefficients) have been based on the charted vertical component 
of the geomagnetic field for 1955.0 (Jensen and Whitaker , 1960; 

Jensen, Murray and Welch , 1960) . Most of the other analyses have 
been based on the more precise charted or observed horizontal components 
of field, for which the fit obtained by Jensen and Whitaker had a 
root-raean-square error estimated to be 1150y (one y = 10 ^ cgs-unit) 
in the charts of the U.S., and 632y in those of the U.S.S.R.; the 
maximum difference in the computed minus observed value of horizontal 
intensity was 6200y. Results based on this analysis have fit, within 
about one per cent, several satellite measurements in the lower 
Van Allen radiation belt (Heppner , et al . , 1960). Coefficients obtained 
by Finch and Leaton have given a similar or somewhat better fit within 
another region (Heppner , et al . , 1960; Finch and Leaton , 1957) . Tabulations 
of field values, of field lines and their conjugate points, and of 
adiabatic invariants applicable to geomagnet ically trapped particles 
have been derived for these coefficients (Jensen, Murray, and Welch , 

I960; Vestine and Sibley , 1960; Ray , et al . , 1962). There is now a 
spherical-harmonic analysis for 1960, based directly on observational 
points ( Jensen and Cain , 1962) . The radiation-belt L-shells of 
Me II wain (1961) use the 1960 values. 
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III . INTEGRALS 


Various books on potential theory have shown how to analyze 

magnetic fields by surface integrals ( Kellogg, 1929). Vestine (1940, 

1941), Taylor (1944), and Benkova (1953) have extended the technique 

to the geomagnetic field of a sphere, while Vacquier , et ah , (1951) 

has done the same for a plane earth. 

A convenient starting point is Green's theorem, which gives a 

magnetic potential V(P) at an internal point P(r,9,\) in terms of 

SV 

surface values of the potential and its normal derivative — (n being 
the outward normal): 

1 I b 1 1 dV \ 

V(P) = - — J I v ) ds » (3) 

4tt S \ 9n r r dn ' 


where S is the surface of the earth. 

If V e is that part of V originating outside the earth, and is 
that part originating inside the earth, then for P external to S, 

1 / d 1 1 dV \ 

v i < p ) - — S v dS > (4 > 

4tt S \ dn r r in ' 

and, for P internal to S, 

1 / d 1 1 dV \ 

V (P) J* V dS , (5) 

e 4tt S \ Sn r r &n / 

where V = V & + V ± . Upon S itself 


V 

e 



1 


- ; 

2tt S 


1 SV 
r dn 


a 1 

(V - U) — - 
dn r 


dS + U, 


( 6 ) 
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where U/4tt is the strength of any uniform double layer on S, a layer 
whose potential is zero outside S, and equal to -U everywhere inside S 


(Ve stine , 1940; Taylor , 1944) 
For a spherical earth, 


1 2tt tt/ 2 

V - V. = — J* J (V + 2a Z) cos i|r di|r dX , 


(7) 


2tt o 


where V - V is taken at the pole of coordinate (a,0,X) and ^ —9/2 
e i 

and Z = dV/dn (Ve stine , 1941) and an analogous expression for Z & - Z ± 
was obtained in correction of Vestine's earlier result (Taylor , 1944). 
Therefore, as in spherical harmonic analysis, an integral method suffices 
to separate an observed surface magnetic field into parts of external 
and internal origin. 

Equation (4) can be transformed with the Green's function of the 
first kind: 


1 a 1 

G(Q,P) , (8) 

R r R' 


where P is the outside point (r,9,\), Q the point (a,0',\') on S , 

R the distance PQ, and R' the distance P'Q (where P' is the image 
point to P in the sphere). From this, we get the well-known Poisson's 
integral : 


V. <r,0,\) 


2 

r 


2 

a 


4Tta 


f (0 ' ,X ') dS 

i “I?— 


(9) 


giving at P outside, in terms of the surface values f(0',\') 
of V. on S. 

l 

Existing spherical-harmonic analyses of the geomagnetic field do 
not show an external contribution that can be detected with any 


- 6 - 


certainty. In this paper, consequently, we will neglect the possible 
contribution of external main-field terms, 

Z is not a potential function, because its direction in space 
is not always the same, but the values X, Y, and Z are known over S 
and are defined in Eq. (2). Therefore, if we take the earth’s center 
0 as origin, we can transform the values into a Cartesian coordinate 
system (x, y, z) with the x axis towards 90°W of Greenwich, the y axis 
toward the Greenwich meridian, and the z axis toward the north pole. 

Thus x * -r sin 0 sin \ , y = r sin 8 cos X, and z = r cos 0, so that 
the transformation results in the potential functions X = - BV/9x, 

Y = - dV/Sy, and Z = - dV/dz. Since X, Y, and Z are potential functions, 
and therefore are not the same as X> Y, and Z, from Poisson’s integral 

r 2 - a 2 X(x',y',z')dS 

X(x,y,z) = J , 

4rra S R J 

r 2 - a 2 Y(x* ,y ' ,z ')dS 

Y(x,y,z) = / 3 , (10) 

4ira S R 

r 2 - a 2 Z(x',y\z')dS 

Z(x,y,z) = J 3 , 

4rra S R 

2 2 2 2 

where x' + y' + z* = a , and X, Y, Z are the surface values of f in 
the x, y, and z directions, respectively. 

The field components on S in the three directions are then 

X ( x '.y' ,z') = X cos e 1 sin X ' - Y cos X ' + Z sin 0 ' sin 1' , 

Y (x ' , y ' , z ' ) = - X cos 0 ' cos X ' - Y sin \ ' - Z sin 0 ' cos X ' , (11) 

ZCx'.y'* 2 ') = X sin 0 - z cos 0 ' , 
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where, as in Eq . (2), 


X 


1 av i av / av 

a 30 1 a sin 9 ' 3X ' \ 3r 



Here, a is the earth's radius, 9' is the colatitude, X' is the 
longitude east of Greenwich, as in Eq. (10). 

By inserting the values of X, Y, an< * Z from charts of the 
geomagnetic field into Eqs . (10) and (11) we may then calculate 
by machine the earth's main field at external points. To calculate 
the field lines and the various adiabatic invariants, we merely 
take the field direction as specified by the field's three orthogonal 
components given by Eq . (10) and use a Runge-Kutta-Gill integration 
scheme (Ve stine and Sibley . 1960) . 

We need know only 3V/3n over S to estimate the potential V at P. 
For a point inside a sphere, P(r,9,X), we use a Green's function of 
the second kind, H(Q,P) . This function gives the potential at P 
relative to a point on the sphere, Q(r', 0',\'), and an image point, 

\ 

— , Q, \ I : 

r ' 1 a 1 2a 2 

H(Q,P)=~+ +“ log — . (12) 

R rR' a a - rr ' cos y + rR' 



The value R is the distance PQ, and the image point P' is inside the 

sphere along the line OP so that the equality of R' and P'Q is 

o 

defined by the condition OP' • OP = a , where a is the radius of 
the sphere (Kellogg . 1929). H(Q,P) can be transformed by known 
methods into a form H'(Q,P) suited to the calculation of V(P) at an 
exterior point P(r,9,X) relative to a point on the sphere Q(r', 0',\'). 
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Thus the equation 

1 

V(P) =— J H'(Q,P) f(Q) dS' , (13) 

4tt S 

provides the magnetic potential at P in terms of the surface values 
dV 

of f(Q) * on S, By differentiating H'(Q,P) with respect to 

r, 0, and and by integrating, we can obtain the field components 

in polar coordinates. But since Eq. (10) seems to offer a more 
complete use of the measured information, we will not use here the 
normal component alone. 

The following sections will describe and discuss for the first 
time new spherical-harmonic analyses of the American and Soviet 
world isomagnetic charts for 1955. We will indicate the quality of 
fit, apply Eq. (10) to the American charts, and compare the computed 
results of field parameters with those obtained by Vestine and Sibley 
(1960) for points in the radiation belts. 
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IV, SPHERICAL -HARMONIC ANALYSES OF ISOMAGNETIC CHARTS 
OF THE U.S, AND U,S,S,R. FOR 1955.0 

Scalings of American and Soviet charts at 5° intervals of 
latitude and longitude have been represented in terms of spherical 
harmonics by one of the authors of this paper (J. L. Carlstedt) . He 
determined the coefficients of Eq. (1) by fitting the values of 
X, Y, and Z, supposing the series to terminate with n = 4, 6, 8, 10, 
or 12. He obtained a weighted least-squares fit for each case by 
using a scheme that Vestine and Lange described in their analyses of 
the isomagnetic charts for 1945 (Vestine , et al . » 1947). 

Table 1 compares the results to with those obtained by Finch 

o 

and Leaton in the British Admiralty charts for 1955. It reveals a 

fairly good agreement in magnitude and in sign, though discrepancies 

-2 0 

are sometimes as great as 10 cgs when results for g^ are quoted 
to four figures. In certain analyses, this means that the computed 
values for the lower Van Allen region may disagree by as much as 
lOOOy (where ly = 10 ^ cgs) . 

We have repeated the analyses, quoting more figures, for the 
American and Soviet charts for 1955. To note the effects dependent 
on the series 1 truncation, we supposed that the series of Eq. (1) 
terminated with P^, P^, pjj, P^, and pjj|. For the American charts 
of X and Y, Tables 2(a) and 2(b) list our values of g^, and h°, 
respectively. Tables 2(c) and 2(d) list our corresponding values 
for Z. Similarly, Tables 3(a), 3(b), 3(c), and 3(d) list our values 
for the Soviet charts. 



- 10 - 


Fanselau and Kautzleben (1958) noted that the terms ud to P; 

6 

4 

(but not to P^) are very roughly the same, whether the series ends 
there or goes on to Agreeing even better, apparently, are the 

results for low-degree coefficients of the series that terminate in 

p 8- p io' and p ir 

Our X, Y, and Z series were synthesized at the points in a 
17 x 36 grid that is composed of 10° spacings in colatitude and 
longitude. Table 4 gives the root-mean-square (rms) errors in fit 
for this grid between the synthesized and charted American and Soviet 
data . 


The quality of fit revealed by the rms errors in the X component 

ranges from about 230y (0.0023 cgs) for the series terminating with 

4 12 

P^ to about llOy for the series terminating in P^ (about 90y for 

the Soviet charts). Interestingly, the series for Z has a larger 

4 

rms departure (320y) when it terminates in P^, and the smallest value 

12 

of all (55y) when it terminates in P^ — a considerable improvement. 
Between the points of a 10° x 10° grid, of course, the present estimates 
of quality will not necessarily apply. In addition, the discrepancies 
for X and Y should theoretically be larger than those for Z . This is 
because we synthesized X and Y from their mean coefficients, not 
from direct analysis, as with Z. In reality, however, the coefficients 
based on Z are probably less accurate than those for X and Y, merely 
because Z has been less accurately measured. Therefore, when Tables 2 
and 3 are used for physical rather than for statistical purposes 
(or for discussion of methodology), it may be appropriate to regard 
the third digit from the right, rounded off, as the last significant 
figure. 
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lt would be interesting to compare the various results obtained 

12 

by using coefficients up to P 12 while computing values in the upper 
atmosphere and beyond. In attempting extrapolations to the earth's 
core, one would probably minimize errors by choosing low-degree terms 
from the series to P^. The analyses contain other defects, some 
not yet mentioned. For instance, the spacing of the 10 x 10 grid 
is necessarily non-uniform in distance. Therefore, though it is 
theoretically possible to determine the coefficients of Eq . (1) so 
that the resulting values are independent of n (the number of terms 
in the series), actual practice may require data at equidistant grid 
points, making difficult the appropriate weighting of data. 

Difficulties also occur in the integral methods mentioned earlier. 
Of principal importance is the distribution of data over the area of 
integration. How this is useful in calculating the geomagnetic field 
will be considered in the next section. But it should be mentioned 
here that non-uniform spacing of sampled chart data could conceivably 
minimize some small effects that might result from the use of uniform 
spacing. Thus, if there are higher-order harmonic components asso- 
ciated with magnetic anomalies, an aliased contribution to lower-degree 
harmonics may arise (Blackman and Tukey . 1959). For analyzing geo- 


physical data, alternative schemes may minimize any possible effects 
(MacDonald, private communication) . 
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V. ANALYSIS BY NUMERICAL INTEGRATION 
TWO EQUAL- AREA GRIDS 

The integral expressions for the field given by Eq. (10) require 

that X, Y, and Z be specified. These integrals can be evaluated 

approximately by using values of X, Y, and Z specified at points on 

a grid covering the sphere, S. Assuming that these field values 

represent average values for surface-area sectors AS, centered on 

k 

the grid points, we can replace the integrals of Eq. (10) with sunmations. 
By selecting grids that specify surface sectors of equal area, we can 
eliminate one source of variation in the integral and use field 
values directly, without weighting. 

Two equal-area grids are described here. The first is formed 
by lines of latitude and longitude, the second by nearly equal 
subdivisions on the faces of a spherical icosahedron. This equal 
spacing of grid points over the surface of the sphere will enable 
similar local features of the surface-field to affect the values of 
X, Y, and Z similarly anywhere on the surface. 

Longitude-Latitude Equal -Area Grid 

A grid of points can be constructed with lines of latitude and 
longitude so that each grid point represents sectors of equal area. 

First, take two planes that are parallel to the earth's equator and 
that are separated by a fixed distance. The area of the earth's surface 
between the planes will be the same, regardless of where the parallel 
planes are located. If we divide the earth's diameter with a series of 
equidistant planes that are parallel to the equator, the spherical 
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ooflBientS between the planes will be equal in area. These spherical 
segments can be further divided into smaller equal areas by drawing 
meridional planes that are spaced at equal polar angles. Figure 1 
shows such a subdivision of the earth's surface. The area sectors 
adjoining the poles are triangular. Each area sector is bounded by 
lines of latitude and longitude. 

To represent this kind of subdivision, we can generate a set 
grid points in a similar fashion. Its advantage is that most 
magnetic field data are similarly scaled at equal intervals of latitude 
and longitude. Its principal disadvantage is the unequal spacing of 
points over the sphere. In a grid with a longitude and latitude 
spacing of 100 km at the equator, for example, the northernmost grid 
points will be 1130 km from the pole, though separated from each other 
by only 18 km in longitude. Giving field data only at these grid points 
would cause unequal representation of the field's surface features due 
to internal sources: the representation of features distributed in 

longitude would obviously surpass that of features distributed in 
latitude. Hence, extrapolations of the field to heights greater 
than 100 km using this grid in combination with Poisson's integral — 
or Eq. (10) — - will reflect more accurately the effects of sources 
that are distributed in longitude than those of sources that are 
distributed in latitude. Nevertheless, extrapolation tests using 
this grid from 70 N to 70 S give results with an accuracy commensurate 
with that of surface data . 

As we shall later use more advanced grid systems, we will here 
use only one calculation as a sample. We take a grid of data points 
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Fig, 1 Method of generating an equal-area grid 
with latitude and longitude lines 
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distributed over the surface of the earth at regular intervals of 
latitude and longitude. By spacing the points closely, we can better 
assess the potential accuracy of the surface integral method. We 
actually use a composite of two latitude-longitude grids. The primary 
grid has points spaced at intervals of 2° in latitude and 4° in longitude. 
The smaller, secondary grid (whose over-all measurements are 4° in 
latitude by 16° in longitude) is located beneath the point at which 
the field is calculated. Its points are spaced at intervals of 

1/4° in latitude and 1/2° in longitude. 

In such a grid, the area sectors assigned to each grid point 
vary in size over the surface of the sphere. We determine the area 
assigned to each grid point by integrating dS = a 2 sin 9 d9 dcp over the 
assigned area, where 6 is the colatitude and cp is the longitude of 
a point on the sphere, and a is the sphere's radius. Values of 
X, Y, and Z for each point were calculated from the coefficients we 
had derived from the American charts of 1955.0; using these values, 
we calculated the field at 50°N, 0°E and at an altitude of 600 km by 
integrating Eq . (10) over the composite latitude- longitude grid. These 
spherical harmonic coefficients give the field at this point as 
0.3654 gauss (with direction cosines of dx/ds = - 0.0723, dy/ds = 0.8937, 
and dz/ds = 0.4464). This value can serve as a reference for 
determining the accuracy of the value given by the surface integral 
(0.3649 gauss with direction cosines of dx/ds - - 0.0723, dy/ds 0-8958, 
and dz/ds = 0.4470). The error in the field's magnitude is about 0.1%. 
The maximum error in the direction cosines is for dy/ds — about 0.25%. 
Obviously, the surface-integral method can yield quite accurate 
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extrapolations of the geomagnetic field. The calculation is quite 
lengthy, however, and requires an interpolation scheme for points at 
which the field is not known. 

Icosahedral Equal-Area Grid 

The foregoing discussion has shown why it is desirable to have 
sectors with equal areas, with a high degree of symmetry about their 
centers, and with shapes that are independent of their location on the 
sphere. From this, J. W. Kern decided that the icosahedron, being 
the highest -degree regular polyhedron, is an excellent model for 
subdividing a sphere . 

If we trace a spherical icosahedron onto the surface of the earth, 
letting one axis correspond to the axis of rotation, we can form an 
equal-area grid by subdividing the faces of the icosahedron with great 
circles. We will regard the grid points so determined as center points 
of the sectors of integration. Figure 2 shows that almost every sector 
will be a hexagon centered on a grid point; the exceptions will be the 
pentagon centered on each vertex of the icosahedron. Figure 2 shows the 
arrangement of the grid points for one face of the icosahedron. Since 
equal areas are to be assigned to each of these grid points, we must 
make small adjustments in the shape of the sectors and distribute them 
over the sphere. Note that the sectors will not, in general, be 
regular spherical polygons. But these irregularities will be quite 
small if we make a fine subdivision, using the technique described below. 

Consider the spherical triangle ABC, shown in Fig. 2. This face 
of the spherical icosahedron is one of the five spherical triangles 
joining in a common vertex at the north pole of rotation. Its sides 
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Spher ical 
icosahedron 



Plane 

icosahed ron 


Fig. 2 Method of generating an equal-area, 
equal-spacing grid based on the spherical icosahedron 
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AB and AC, lie in meridional planes and can be subdivided into m 
equal arcs. (In the figure, m — 16.) This subdivision determines 
points on AB and AC that can be identified by their colatitude 
measured from the pole. If we pass great circles through the points 
on AB and AC that have the same colatitude, we will determine m 
great-circle arcs inside the ABC face of the spherical icosahedron. 

MN in Fig. 2 is one such arc. If we now subdivide each arc once for 
every arc lying between it and the pole, we generate a regular grid. 

This grid can be transferred to all twenty of the spherical triangles 
that make up the spherical icosahedron. For example, the grid for 
triangle ABC in Fig. 2 may be transferred to triangle BCD by rotating 
it 180 about a radius through E, the midpoint of BC . By joining two 
spherical triangles (such as ABC and BCD in Fig. 2) we form an equilateral, 
spherical quadrilateral (ABDC) . To cover the sphere, we need only 
repeat the grid pattern of ABDC five times around each of the poles. 

Thus, identification of this single pattern will enable us to prescribe 
the grid for an entire sphere. 

This technique can be formalized in terms of a set of vector 
operations, in which linear transformations begin by generating the 
grid within the spherical quadrilateral, and then repeat the grid 
pattern about each pole. It can be made so general that an electronic 
computer can readily generate the coordinates for a grid of any 
frequency, m. 

NUMERICAL FORMULAE FOR SURFACE INTEGRATION 

After a grid has been generated, and after the Cartesian components, 

X, Y, and Z, have been calculated, the integral expression for a given 
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component at any point of space can be replaced by this summation: 



Here, X(x,y,z) is the desired x-component of the field in space 
(similar expressions can be written for Y and Z) ; k is the index of 
summation; X k (x^, y^, z^) is the field value for AS fc (the k-th surface- 
area sector) that is specified by the components of the radius vector 
to the center point (x^, y' and z£) ; and is the distance from the 
point in space to the center of the k-th sector . 
is given by 


R 2 = r 2 + a 2 - 2(xx^ + yy^ + zz^) , 

where r is the radius vector from the center of the sphere to the point 
in space, a is the radius of the sphere, and (x,y,z) and (x^, y^, z^) 
are the Cartesian coordinates of the point in space and of the k-th 
sector's center. 

Both the longitude-latitude and the icosahedron grids assigned 

equal areas to the centers of their surface-area sectors. Therefore, 

2 2 

each sector is AS = 4na /N, where 4-na is the area of the sphere, and 

N is the total number of grid points. If we substitute this expression 

for AS, into the summation for X given above, we see that 
k 

a (r 2 - a 2 ) N » 

X(x,y,z) = £ -r , (14) 

N k=l r£ 
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while 

a(r 2 - a 2 ) N Y 

Y(x,y,z) = £ — , 

N k=l RjJ 

and 

a (r 2 - a 2 ) N Z 

Z(x,y,z) = £ . 

N k=l RjJ 

These are the approximate expressions that will be used in the surface- 
integration method of calculating the field in space from surface values 
over the earth. For example, the total number of grid points over the 

n 

surface of a spherical icosahedron is 10m + 2, where m is the number 

of subdivisions on the icosahedron's edges. 



to reduce the computation's complexity. For a point near the surface 
of the earth, R becomes small for some terms in the series. As R 
approaches zero, the integrand approaches infinity. This challenges 
the range of applicability for Eq. (15) and similar approximations. 

If X^ is a good representation of X over AS. , we may write 

sC 
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XdS 


dS 


Jj~‘‘ J ? 


AS k R ““k 


and then estimate 


S5 k R 

the appropriate variable weight by this equation: 


dS 

« - J — 

AS k R 


One such value is, of course, 


W = 


4rra 

K 


(16) 


An alternative would 
A0 (Fig. 3): 


be the integral for the cap segment of angle 


2tt A9 


sin 0 d8 dcp 


W J 0*^ (a 2 + r 2 - 2ar cos 9) 3 ^ 2 


2tt A© 

- I J 


sin 0 d0 dcp 


0 0 V 2 + z 2 + a 2 -2(x a cos cp sin 0 + z Q a cos 0)j 

Loo ° 


3/2 


For small A0, the approximations 


0 * 

sin 9 ~ 8> cos ® “ 1 

2 


after making a - 1, lead to the integral 


2tt A0 

J S 


0 d0 dcp 


j w p „ 3/2 

0 o |x 2 + z 2 + 1 - 2z - 2x 0 cos cp + z o 9 J 
L O o O 


( 17 ) 
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Z 



Fig. 3 Coordinate system for surface integration 
over a single area sector 
(circular cap on sphere) 


-23- 


The integration of Eq . (17), with respect to 8, yields 
1 2rr |k 2 (K 2 -2x q A8 cos cp + z p A8) 2 - (K 2 -x q A6 cos cp)^j <kp 

_ J — — ~ 

K 0 ( K 2 Z _ x 2 cog 2 C p)( K 2 -2x A0 cos cp + z A0) 

N o o ° 

where 

2 2 . 2 
K = 1 - 2z o + X o Z o 

Table 5 exhibits the relative merits of Eqs . (16) and (18) for 

the case x =0. In it, A0 = .04, which is quite close to the weighting 
o 

appropriate to N = 2562. 

For x t 0, Eq. (18) must be evaluated by numerical quadrature. 

O 

In subsequent numerical results, an eight-point Gaussian quadrature 
was used over (0 ,tt) . Figure 4 exhibits the behavior of one-half of 

2 2 h 

the integrand, both for values of x q and z q where (x q + Z Q ) “ 1,l » 
and for positions ranging from directly over the center of AS k out 

to a distance of one sphere radius . 

A trapezoidal integration of the curve for D - 1 (Fig. 4) yields 

a value of .00504. 

Comparing this result with corresponding entries in Table 5, we 

see that for a distance of one sphere radius from the center of AS^, 

Eqs. (16) and (18) give virtually the same result regardless of 

orientation. This suggests that a combination of Eqs. (16) and (18) 

satisfactorily estimates the surface integral. Equation (16) should 

be used for those AS k farther than one sphere radius from the point 

in space, and Eq . (18) for those closer than one radius. 

dS 

Table 6 estimates J -r for a grid corresponding to N - 2562. 

S R j 

The Center column lists positions along a radius vector through the 
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center of a particular AS^; the Off Center column lists positions 

along a radius vector through a point midway between two AS^> To 

account for overlapping AS^ and to achieve the distribution of errors 

listed under Center, A0 in Eq . (18) and the coefficient in Eq . (15) 

2 

were modified as follows: A9 = (.0675); 

4 „ VN~ 

the coefficient = — (.9978). 

N 

In general, one-tenth of a sphere radius is the minimum distance 
from the sphere at which an error of less than 17 o can be maintained. 
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VI. APPLICATION TO GEOMAGNETIC FIELD-LINE CALCULATION 

To avoid reading the magnetic charts at the previously described 
grid points and to get results comparable to the spherical harmonic 
analyses, we obtained the surface values of X, Y, and Z from the 
48-coefficient expansion of the American 1955 charts. These data 
allow one to compare directly the field lines obtained by the integral 
method with those lines obtained by Vestine and Sibley (1960) * 

Each integration was by a fourth-order Runge-Kutta-Gill scheme, 
moving by steps of 0.1 earth radius ( Gill , 1951). 

Table 7 and Fig. 5 indicate the results of generating a field line 
starting at 0.1 of a radius above the point at 40°N, 90°E and ending 
at an altitude of about 0.1 of a radius. 

Figure 6 shows a 50°N, 0°E field line generated by both the spherical- 
harmonic and the surface-integral methods. For each, the integration 
step is 0.25 earth radius. Note that the difference between the 
calculated positions in the equatorial plane is greater than in Fig. 5 
where we have a smaller step-size. Near the surface, the Runge-Kutta- 
Gill integration of Eq. (16) yields field-line positions quite close 
to those obtained by the spherical-harmonic analysis. Differences in 
conjugate points at constant altitude are less than 0.01 earth radius 
(about 65 km). To determine mirror points (the point along the line 
where the field is the same as at the starting point), the integration 
of Eq. (16) is less useful. The calculated magnitude of the field can 
have errors on the order of 3% at altitudes of less than 0.1 earth 
radius. Such an error would allow the altitude of mirror points 
estimated with Eq. (16) to be in error by as much as 0.01 earth radius. 
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The field lines calculated with either Eq . (16) or with the 
two-part integral correspond quite closely to those calculated by 
spherical-harmonic methods. This correspondence can be further improved 
by using shorter integration steps in the surface-integral methods. 

This would require more computing time than in the calculations of 
Figs. 5 and 6. But because surface-integral methods already increase 
more than tenfold the time required by the spherical-harmonic methods, 
the penalty for using shorter steps is proportionally small. 

Table 7 lists the coordinates for three lines. The first was 
derived by Vestine and Sibley (1960); the second by the mixed- integra 1 
method discussed above; and the third by the integral method that uses 
Eq . (15) throughout. For each step along these three lines, we have 
computed the magnetic-field values (X,Y,Z) and listed them in Table 8. 

The three methods produce quite comparable results for shorter 
field lines. Figure 6 and Tables 7 and 8 show that the two integral 
methods produce almost identical results that differ from the spherical 
harmonic ones by about 100 km. The consistency of the integral results 
can be explained in part by the fact that, although the values of X, 

Y, and Z in Table 8 differ, they are almost proportional. The constant 
of proportionality is the square of the ratio of the factors noted 
in Table 6. This is not entirely unexpected. Because the direction 
of the line of force is determined by the direction cosines associated 
with the X, Y, and Z at each point, the constant of proportionality 
disappears . 


In general, as a point moves away from the surface of the earth, 

Eq. (15) provides an increasingly good approximation of the field. Near 
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the surface, however, only a few values tend to be heavily weighted, 

3 

because of the 1/R argument in Eq. (15). This weight disappears in 
the computation of the direction cosines. A clue to this behavior may 
be found by rewriting Eq. (14): 



3 3 3 3 3 

where ^ ” R x R 2 ' ' * ^k-l ^k+l *'* T ^ e direction cosine defined by 

X becomes 





precisely the direction cosine at the surface. 
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VII. COMPARISON OF CONJUGATE POINTS CALCULATED FROM 
VARIOUS SPHERICAL-HARMONIC ANALYSES 

DeWitt's study of IGY all-sky camera data (DeWitt, 1962) has 
shown that auroras occur simultaneously in regions that are connected 
along magnetic field lines (that is, geomagnetically conjugate). 

Vestine and Sibley (1960) computed the conjugate points of a number 
of locations in the auroral zones. They calculated that the IGY all-sky 
camera station at Farewell, Alaska, and the station at Campbell Island 
in the Southern Pacific Ocean are nearly conjugate, and that the same 
is true of Kotzebue, Alaska, and Macquarie Island. DeWitt's observations 
indicate that Vestine and Sibley may have been in error by as little 
as 20 km. 

We have repeated their calculations for Campbell Island, using 
spherical-harmonic coefficients to order 6. Table 1 lists these 
coefficients we derived from the British (Finch and Leaton , 1957), American, 
and Soviet isomagnetic charts for 1955.0. We also used coefficients 
derived from the American charts for I960 (Jensen and Cain , 1962). 

We have extrapolated the field line from Campbell Island (52°32'S, 

168°59'E) to the northern hemisphere by the method of Vestine and Sibley 
(1960). The surface positions of points conjugate with Campbell Island, 
calculated with the indicated coefficients are: (1) British, 61°59'N, 

154°39'W; (2) American, 1955, 61°45'N, 156°46'W; (3) Soviet, 61°40'N, 
156°23'W; (4) American, 1960, 61°37'N, 155°05'W. The mean latitude 
of the four conjugate points is 61°45'N, with a mean deviation for the four 
positions of ± 7' (about ± 13 km). The mean longitude of the four conjugate 
points is 155°43'W, with a mean deviation of ± 52' (about ± 44 km at the 
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mean latitude). The probable error in this procedure has been 
previously estimated to be about 8 km. Apparently, uncertainties in 
the mathematical representations of the surface field can cause 
significant errors in the conjugate points. For the mean latitude 
and longitude of the four conjugate points of Campbell Island, the 
probable errors are ±4' = ± 7 km and ± 26' = ± 22 km. 

DeWitt's observations (1962) suggest an error limit of about this 
order for the Campbell Island calculations of Vestine and Sibley (I960). 
Obviously, no absolute criterion exists for preferring any one of the 
above determinations. Note that the largest probable error is in the 
conjugate point's longitude. Owing to the general east-west elongation 
of auroral structures, the longitude may also be the most difficult to 
check by auroral observations. 


-33- 


VIII, MAGNETIC MOMENT AND DIPOLE AXIS SINCE 1833 

Spherical-harmonic analyses of British, American, and Soviet 

isomagnetic charts for 1955.0 give coefficients that agree rather well 

with each other. This is not surprising; they are independently derived 

from practically the same sets of measurements. Truncating the series, 

however, does affect the results. We see in Fig. 7(A) the positions 

of the geomagnetic north pole since 1829, as placed by series to 4, 

to 6, and to 12. It will be seen that pole positions for 1955 from 
4 

our analyses to agree well with those of Gauss (1839) for 1838. Also 
agreeing well with each other are those positions derived from series 
to P^, including the last analysis (Jensen and Cain . 1962). 

Figure 7(B) shows that the earth's dipole moment still continues 
its rather uniform decrease with time. Points for Soviet charts are 
based on our analysis presented here, as well as on a recent analysis 
by Adam , et al . . (1962) . 




Magnetic moment, M (cgs) 



Yeor 


Fig. 7 (A) Geomagnetic north pole, from various analyses, 1829 — 1960 

25 

(B) Moment, M, in units of 10 cgs 
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IX. MAGNETIC FIELD IN THE EARTH'S MANTLE 

As we have seen, only a very small part of the magnetic field at 
the earth’s surface can be ascribed to external sources. The main field 
is, in fact, thought to originate in the earth's core, and to be 
maintained by motions of electrically-conducting fluid material in the 
core. Extrapolation of the surface field to the boundary of the core 
is therefore of great interest to geophysicists. However, the sources 
producing the surface field are not entirely confined to the core. The 
spherical-harmonic representation of the field given by Eq . (1) may 
contain contributions from sources in the earth's mantle and crust. 

For illustration, we have calculated the nondipole portion of the main 
field at a depth of 1500 km (about halfway down to the core), using 
the American coefficients of 1955.0. 

Figure 8 shows the results in terms of the downward (Z) and 
horizontal (H) components of the field. To exclude contributions from 
localized sources in the outer mantle and crust, it uses only terms 
up to m=n=6. If such contributions were present, the strong dependence 
of higher-order terms on r in Eq . (1) would produce relatively large 
errors in the field calculated for r < a. This kind of extrapolation 
seems to present no numerical difficulties, and if the possible sources 
in the crust and mantle are reflected only in coefficients of degree 
and order higher than 6, the extrapolation can in principle be 
extended to the boundary of the core. 
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Fig. 8 Nondipole main-field components for 1955.0 at a depth of 1500 km 
are for downward (£) component in gauss; horizontal component (H) 

is given by arrows 


Contours 
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X . CONCLUSIONS 


Spherical-harmonic or integration analysis of the geomagnetic 

field's potential offers a way to extrapolate and interpolate the field 

over the earth's surface and nearby space. Spherical-harmonic analyses 

of different sets of isomagnetic charts yield sets of coefficients that 

are nearly the same to about with differences commensurate with the 

o 

probable errors in the data. To represent complex surface features 
(such as anomalies), we must use a spherical-harmonic series with many 
higher-order terms. Theoretically, integral methods avoid this difficulty. 
But if accuracy comparable to spherical-harmonic analysis is required, 
similar difficulties are met in supplying field values for the large 
number of data points that are necessary for an integration. 

Geomagnetic field lines can be extrapolated into space either by 
spherical-harmonic analyses or by surface integrals. To this 
extrapolation, Poisson's integral is applied. The key is recognizing 
that each component of the geomagnetic field, specified over the earth's 
surface in a Cartesian-coordinate frame rotating with the earth, can 
be treated as a scalar potential function. Then, by applying Poisson's 
integral separately to each component, one can calculate the corresponding 
components at any point in space external to the sphere. We doubt that 
such application to a sphere has been previously noted. 

The calculations based on ‘spherical harmonics take less computation 
time, because their coefficients embody the surface-field data in a 
convenient analytic form. Surface-integral calculations require a 
relatively complete specification of the surface field, particularly 
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those parts near the earth's surface. Each integration over the surface 

requires the handling of the total surface-field data. It involves an 
3 

argument of 1/R , where R is the distance from a point at which the 
integral is evaluated to a data point on the surface of the sphere. 

The calculation of this for each data point greatly extends the time 
required for evaluating the surface integral. Accuracies comparable 
to those with spherical-harmonic methods require increases in time by 
factors of between ten and one hundred. 

However, since the integral method involves only the process of 
summation, it may be less limited in Accuracy. It also offers 
conceptual advantages convenient to those familiar with potential 
theory. Further, the calculations presented here indicate that 
comparable field-line positions are obtained by both methods. 

There is no evidence of change in the position of the geomagnetic 
north pole from 1833 to I960. The apparent motions are in some 
instances probably due to truncated series in spherical -harmonic 
representations of the field. The dipole moment of the earth has 
decreased at a uniform rate from 1829 to I960. 

Finally, there appear to be no serious obstacles to an inward 
extrapolation of the geomagnetic field to the earth’s core. Sources 
of the geomagnetic field external to surface of the earth are negligi- 
ble, and those due to sources in the mantle and crust are likely to 
contribute only to harmonics of order and degree higher than 6. Thus 
a truncated spherical-harmonic series should give an adequate represen- 


tation of the field external to the core. 
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APPENDIX 


Table 

1. Spherical-Harmonic Analyses of the Geomagnetic Field 

for 1955, Based on British, American, and Soviet 
World Charts. Gaussian (Schmidt) Coefficients, 

§n’ n> in units cgs 

2. Spherical-Harmonic Analysis of the Geomagnetic Field 

for 1955, U.S. World Charts. For Series Termin- 
ating with Indicated Values in P^ 

a. Mean Values of X and Y, Coefficients g m . . . 

— — > °n 

b. Mean Values of X and Y, Coefficients h m . . . 

c. Mean Values of Z, Coefficients g 

— 5 °n 

d. Mean Values of Z, Coefficients h m 

3. Spherical-Harmonic Analysis of the Geomagnetic Field 

for 1955, U.S.S.R World Charts. For Series 
Terminating with Indicated Values in P^ 

a. Mean Values of X and Y. Coefficients g m . . . 

b. Mean Values of X and Y, Coefficients h m . . . 

c. Mean Values of Z y Coefficients g n 

d. Mean Values of Z, Coefficients h m 

— 5 n 

-2 

4. Weighted rms Error E in Units of 10 cgs in the Fit 

of Tabular Values of 10 Intervals of Colatitude 
and Longitude; U.S. and J.S.S.R. Charts for 1955 
for Spherical-Harmonic Series Terminating in 
Degree 4, 6, 8, 10, and 12 for (a) U.S. Charts 
and (b) U S.S.R. Charts 


5. 

6 . 

7. 


Values of Eqs , (16) and (18) for A9 = 0.4 and 


x = 0 
o 


Values of J — for Various Heights Above the Sphere 


S R 


Position Coordinates Estimated by Three Methods. 
Field Line Starts at 40°N, 90°E; h * 0.1 Earth 
Radius 


8. Components of Geomagnetic Field Estimated by Three 
Methods. Field Line Starts at 40°N, 90°E; 
h 53 0.1 Earth Radius 
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Table 1 

Spherical-Harmonic Analyses of the Geomagnetic Field for 1955, 
Based on British, American, and Soviet World Charts. 
Gaussian (Schmidt) Coefficients, g™, h§, 
in units of 10”^ cgs 



m 

g n 

h* 

n 

n,m 

British 

U.S. 

U.S.S.R. 

British 

U.S. 

U.S.S.I 

I - .... 

1,0 

-3055 

-3055 

-3051 




2,0 

- 152 

- 147 

- 141 




3,0 

118 

117 

113 




4,0 

95 

87 

97 




5,0 

- 27 

- 24 

- 33 




6,0 

10 

2 

7 




1,1 

- 227 

- 210 

- 202 

590 

585 

584 

2,1 

303 

307 

299 

-190 

-185 

-187 

3,1 

- 191 

- 170 

- 174 

- 45 

- 59 

- 56 

4,1 

80 

65 

78 

15 

18 

11 

5,1 

32 

40 

33 

2 

10 

10 

6,1 

5 

12 

8 

- 2 

- 6 

- 7 

2,2 

158 

145 

168 

24 

49 

38 

3,2 

126 

127 

125 

29 

30 I 

26 

4,2 

58 

47 

57 

- 31 

- 24 

- 31 

5,2 

20 

21 

16 

10 

10 

12 

6,2 

2 

3 

3 

11 

16 

14 

3,3 

91 

86 

81 

- 9 

- 3 

- 8 

4,3 

- 38 

- 44 

- 36 

- 4 

- 7 

- 5 

5,3 

- 4 

5 

9 

- 5 

- 1 

- 4 

6,3 

- 24 

- 26 

- 26 

0 

0 

- 1 

4,4 

31 

29 

33 

- 17 

- 13 

- 13 

5,4 

- 15 

- 15 

- 14 

- 14 

- 18 

- 14 

6,4 

- 3 

7 

- 3 

- 1 

0 

0 

5,5 

7 

- 4 

6 

9 

8 

1 

6,5 

0 

3 

1 

- 3 

- 4 

- 1 

6,6 

- 11 

- 10 

9 

- 1 

- 4 

- 3 















is of the Geomagnetic Field for 1955, 
of X and Y , for series terminating 


ients, (Units of 10 cgs) 


*TJ 

00 00 

10 

10 

p!2 

12 

0 

0 

0 

-305191 

-305179 

-305189 

- 2227T 

-21996 

-22253 

-14735 

-14740 

-1471*0 

30705 

30727 

30734 

13957 

13750 

13578 

12179 

12225 

12235 

-17411 

-17287 

-17386 

12722 

12723 

12702 

9155 

9129 

9327 

8677 

8620 

86l6 

6541 

6578 

6592 

4484 

4385 

4299 

-4499 

-4531 

-4544 

3085 

3342 

3485 

-1677 

-1581 

-1519 

3683 

3777 

3701 

2081 

2083 

2048 

-176 

-191 

-85 

-1716 

-1759 

-1806 

-692 

-644 

-615 

201 

83 

65 

1322 

1372 

1394 

-409 

-479 

-546 

-2723 

- 2766 

-2786 

-577 

-439 

-359 

218 

218 

212 

-1241 

-1380 

-1470 

971 

1118 

1248 

-634 

-554 

-620 

79 

81 

32 

152 

141 

217 

-328 

-380 

-441 

-407 

-379 

-360 

196 

136 

118 

310 

378 

307 

-31 

-215 

-249 

10 

65 

96 

-103 

-160 

-216 

-268 

-325 

-352 

79 

175 

232 
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®ible 2(a)eont'd 



— 

Gaussian (Schmidt) Coefficients, g (Units of 10 cgs) 

n, m 

< 

o 

*8 

P 10 

10 

P 12 

12 

8, 5 



2 

4 

-5 

8, 6 



-157 

-240 

-293 

8 f 7 



116 

59 

28 

8, 8 



163 

216 

204 

9 , 0 




197 

392 

9 , 1 




-168 

-228 

9, 2 




-29 

-94 

9 > 3 




-44 

17 

9, 4 




-36 

-108 

9, 5 




109 

126 

9, 6 




-44 

-64 

9, 7 




95 

51 

9, 8 




30 

2 

9, 9 




33 

42 

10, 0 




-241 

-291 

10, 1 




-96 

-60 

10, 2 




-97 

-145 

10, 3 




-159 

-191 

10, 4 




88 

131 

10, 5 




19 

7 

10, 6 




-98 

-138 

10, 7 




-96 

-135 

10, 8 




11 

3 

10, 9 




51 

66 

10, 10 




12 

7 

11, o 





255 

U, 1 





117 

11, 2 





-255 

11, 3 





136 

11, 4 





-50 

11, 5 





93 

11, 6 





-2 

11, 7 





-34 

11, 8 





-41 

11, 9 





2 

11, 10 





-31 

n, n 





1 

12, 0 





-62 

12, 1 





7 

12, 2 





85 

12, 3 





-74 

12, 4 





4 

12, 5 





-80 

12, 6 





-73 

12, 7 





-49 

12, 8 





23 

12, 9 





-2 

12, 10 





-34 

12, 11 





-4 

12, 12 


] 

. 

- 

0 
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Table 2(b) Spherical -Harmonic Analysis of the 'Vomagnetic Field for 1955, 

U. S. World Charts. Mean values of X and Y , for series terminating 

with Indicated values in I ,:n . 

n 
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Ssble 2(b) cont ' d 



— — ■ ■ ■■ — ■■ . ■ - — - r ■■ ■ 11 — 

Gaussian (Schmidt) Coefficients, h Q (Units of 10” cgs) 


_4 

_6 

_8 

„10 

_12 

n, m 

P 4 

P 6 

? 8 

P 

10 

P 

12 

8, 5 



118 

68 

72 

8 1 6 



91 

11 

-29 

8, T 



42 

92 

100 

8, 8 



-282 

-279 

-330 

9, 0 




0 

0 

9, 1 




-255 

-230 

9, 2 




251 

336 

9, 3 




174 

148 

9, U 




9 

27 

9, 5 




128 

194 

9, 6 




163 

166 

9, 7 




120 

27 

9, 8 




-30 

-24 

9, 9 




-65 

-42 

10, 0 




0 

0 

10, 1 




231 

191 

10, 2 




48 

43 

10, 3 




-134 

-176 

10, 4 




-17 

-34 

10, 5 




-104 

-101 

10, 6 




27 

-3 

10, 7 




64 

73 

10, 8 




16 

-16 

10, 9 




19 

25 

10, 10 




-6 

-9 

11, 0 





0 

11, 1 





-155 

n, 2 





49 

11, 3 





103 

11, 4 





26 

n, 5 





75 

11, 6 





-12 

11, 7 





-101 

11, 8 





35 

11, 9 





10 

11, 10 





37 

n, n 





0 

12, 0 





0 

12, 1 





27 

12, 2 





130 

12, 3 





-203 

12, 4 





60 

12, 5 





-37 

12, 6 





-48 

12, 7 





-9 

12, 8 





-40 

12, 9 





2 

12, 10 





21 

12, 11 





-2 

12, 12 




_ 

0 
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Table 2(c) Spherical-Harmonic Analysis of the Geomagnetic Field for 1955 » 
U. S. World Charts. Mean values of 1 , for series terminating 

with Indicated values in P™. 


Dj m 

Gaussian (Schmidt) Coefficients, g™ (Units of 10”^ cga) 

A 

A 

4 


P 12 

12 

0, 0 

3 b 

321 

245 

197 

134 

1, 0 

-302261 

-304385 

-304137 

-304074 

-303936 

1, 1 

-20508 

-19510 

-19530 

-19542 

-19552 

2, 0 

-15723 

-15231 

-15421 

-15528 

- 15616 

2, 1 

28576 

28769 

28915 

28915 

28916 

2, 2 

16545 

16581 

16660 

16668 

16660 

3, 0 

13593 

11144 

11506 

11629 

11800 

3, 1 

-21123 

-18413 

-18463 

-18517 

-18506 

3, 2 

10791 

11485 

11319 

11252 

11250 

3, 3 

10166 

10025 

10016 

9988 

9974 

k , 0 

8723 

9315 

9026 

8862 

8764 

4, 1 

7158 

7492 

7773 

7776 

7752 

4, 2 

5459 

5609 

5859 

5893 

5895 

4, 3 

-319k 

-3735 

-3697 

-3696 

-3694 

k , 4 

3719 

3651 

3655 

3659 

3660 

5, 0 


-2663 

-2192 

-2017 

-1773 

5, 1 


3969 

3893 

3773 

3844 

5, 2 


1875 

l48o 

1337 

1318 

5, 3 


-752 

-811 

-923 

- 978 

5, 4 


-1098 

-1120 

-1116 

-1122 

5, 5 


-220 

-252 

-254 

-256 

6 , 0 


687 

323 

96 

-32 

6 , l 


445 

890 

895 

791 

6 , 2 


300 

750 

827 

862 

6, 3 


-1952 

-1827 

-1826 

-1813 

6, 4 


-450 

-407 

-390 

-385 

6, 5 


109 

127 

143 

154 

6 , 6 


-558 

-552 

-552 

-552 

7, 0 



538 

757 

1078 

7, 1 



-103 

-291 

-140 

7, 2 



-647 

-896 

-939 

7, 3 



-157 

-396 

-524 

7, 4 



-100 

-82 

-96 

7, 5 



-247 

-275 

-287 

7, 6 



-73 

-95 

-94 

7, 7 



521 

529 

528 

8, 0 



-4l8 

-695 

-864 

8, 1 



598 

607 

407 

8, 2 



6 99 

829 

919 

e, 3 



254 

256 

285 

8, 4 



146 

189 

201 
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Ifcble 2(c) cont'd 


' " ' ' “ X" — 1 

Gaussian (Schmidt) Coefficients, g“ (Units of 10“ cgs) 


n, m 


A 

A 

P 10 

10 

12 

12 

8, 5 



93 

206 

234 

8, 6 



6o 

45 

55 

8, 7 



-116 

-125 

-124 

8, 8 



137 

137 

137 

9, 0 




247 

643 

9, 1 




-256 

-24 

9, 2 




-372 

-447 

9, 3 




-432 

-662 

9 , 4 




44 

14 

9, 5 




-118 

-153 

9» 6 




-131 

-124 

9, 7 




86 

62 

9, 8 




-100 

-107 

9> 9 




-29 

-26 

10, 0 




-316 

-524 

10, 1 




12 

-289 

10, 2 




194 

348 

10, 3 




2 

55 

10, 4 




92 

n6 

10, 5 




343 

414 

10, 6 




-80 

-4o 

10, 7 




-6l 

-14 

10, 8 




-7 

4 

10, 9 




-no 

-n2 

10, 10 




-78 

-72 

n, o 





449 

n, i 





305 

11, 2 





-109 

n, 3 





-374 

11, 4 





-56 

11, 5 





-90 

11, 6 





23 

11, 7 





-137 

11, 8 





-50 

11, 9 





34 

11, 10 





-48 

n, n 





-32 

12, 0 





-248 

12, 1 





-376 

12, 2 





232 

12, 3 





82 

12, 4 





43 

12, 5 





152 

12, 6 





n8 

12, 7 





200 

12, 8 





73 

12, 9 





-14 

12, 10 





81 

12, 11 





10 

12, 12 





1 
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Tfcble 2(d) Spherical -Harmonic Analysis of the Geomagnetic Field for 1955, 

U- S. World Charts. Mean values of Z , for series terminating 

with indicated values in P®. 

n 


n, m 

Gaussian ( 

Schmidt) Coefficients, h'^ (Units of 10”^ cgs) 


p 6 

4 

| F 10 

10 

*12 

0, 0 

0 

0 

0 

0 

0 

1, 0 

0 

0 

0 

0 

0 

1, 1 

56779 

56787 

56736 

56739 

56728 

2, 0 

0 

0 

0 

0 

0 

2, 1 

-18000 

-18067 

-18115 

-18117 

-18117 

2, 2 

4418 

4575 

4569 

4572 

4564 

3, 0 

0 

0 

0 

0 

0 

3, l 

1 -4336 

-4315 

-4440 

-4424 

-4412 

3, 2 

1 1749 

2042 

2041 

2073 

2073 

3 > 3 

i -544 

-551 

-525 

-510 

-520 

k , 0 

! ° 

0 

0 

0 

0 

4, 1 

i 1036 

922 

829 

8o4 

802 

h , 2 

1 -2889 

-2230 

-2248 

-2236 

-2234 

4, 3 

i -1324 

! -1375 

-1428 

-1429 

-1426 

1+, 4 

-1876 

-1942 

-1948 

-1949 

-1952 

J f 0 


! 0 

0 

0 

0 

2, 1 


32 

-157 

-123 

-43 

5.- 2 


792 

790 

858 

857 

5, 3 


-36 

124 

187 

148 

4 


-1014 

-923 

-919 

-913 

5, 5 

6 , 0 


580 

0 

562 

0 

565 

0 

562 

0 

6, 1 


-153 

-299 

-353 

-362 

6 , 2 


1320 

1288 

1314 

1353 

6, 3 


-182 

-354 

-357 

-340 

6, 4 


-430 

-489 

-495 

-504 

6, 5 


114 

186 

177 

179 

6, 6 


-451 

-456 

-455 

-455 

7, 0 



0 

0 

0 

7, 1 



-256 

-202 

-32 

7, 2 



-3 

114 

112 

7, 3 



426 

559 

469 

7> 4 



4o8 

426 

442 

7, 5 



-139 

-104 

-124 

7, 6 



-69 

-57 

-60 

7, 7 

8 , 0 



-344 

0 

-338 

0 

-337 

0 

8, l 



-197 

-282 

-301 

8, 2 



-50 

-5 

91 

3 



-348 

-353 

-317 

8 , 4 



-200 

-217 

-237 
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Table 3(a) Spherical -Harmonic Analysis of the Geomagnetic Field for 1955, 

U.S.S.R. World Charts. Mean values of X and Y , for series terminating 

vith indicated values in P®. 

n 



Gaussian (;*: 

'chraidt) Coefficients, g™ (Units of 10 ^ cgs) 


4 

_6 

_8 

1 r -)10 

12 


P 4 

P 6 

P 8 

P 

10 

P 12 

n, m 






0, 0 

0 

0 

0 

0 

0 

1, 0 

-304078 

-305135 

-304975 

-304965 

-304973 

1- 1 

-28474 

-20227 

-21230 

-20832 

-21195 

2, 0 

-14495 

-14138 

-14150 

-14154 

-14154 

2, 1 

29756 

29877 

29909 

29926 

29935 

2, 2 

15913 

16781 

16814 

16329 

15898 

3, 0 

13593 

11250 

11571 

11609 

11618 

3, 1 

-20413 

-17440 

-17793 

-17641 

-17776 

3, 2 

11843 

12450 

12425 

12472 

12485 

3, 3 

9222 

8093 

8554 

8207 

8236 

4, 6 

9127 

9690 

9668 

9629 

9624 

4, 1 

7626 

7818 

7875 

7904 

7916 

4, 2 

5326 

5723 

5726 

5497 

5298 

4, 3 

-2499 

-3574 

-3625 

-3641 

-3656 

4, 4 

3427 

3248 

3404 

3395 

3641 

5 , 0 


-3255 

-2795 

-2715 

-2664 

5i 1 


3259 

3044 

3148 

3052 

5, 2 


1626 

1596 

1666 

1688 

5, 3 


- 899 

-666 

-845 

-829 

5j U 


-1423 

-1621 

-1601 

-1607 

5, 5 


-581 

-882 

-674 

-872 

6. 0 


725 

691 

609 

592 

6 , 1 


837 

907 

947 

966 

6, 2 


342 

345 

185 

46 

6. 3 


-2623 

-2688 

-2710 

-2731 

6, 4 


-312 

-225 

-228 

-90 

6, 5 


139 

176 

135 

117 

6, 6 


-864 

-1114 

-1050 

-1109 

7. 0 



609 

732 

839 

7# 1 



-773 

-687 

-768 

7# 2 



296 

389 

421 

7, 3 



134 

10 

21 

7* H 



-413 

-385 

-393 

7. 5 



- 297 

-176 

-287 

7, 6 



297 

258 

262 

7* 7 



206 

153 

210 

8. 0 



-45 

-172 

-205 

8. 1 



-18 

29 

56 

8. 2 



-553 

-682 

-795 

8. 3 



-57 

-84 

-112 

8, 4 



126 

125 

224 
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T&ble 3(a) cont’d 


n, m 

Gaussian (Schmidt) Coefficients, g“ (Units of 10” b cgs) 

K 

4 

08 

*8 

P 10 

10 


0, 5 



122 

73 

50 

8 , 6 



-212 

-174 

-209 

0 , 7 



102 

32 

22 

8, 8 



213 

163 

252 

9 , 0 




164 

325 

9 , 1 




88 

14 

9 , 2 




277 

320 

9 , 3 




-156 

-147 

9 , 4 




89 

79 

9 , 5 




167 

89 

9 , 6 




-79 

-74 

9 , 7 




-14 

21 

9 , 0 




13 

34 

9 , 9 




-6 

-7 

10, 0 




-166 

-215 

10, 1 




7 

46 

10, 2 




-354 

-454 

10, 3 




-23 

-52 

10, 4 




25 

103 

io, 5 




-29 

-57 

10, 6 




05 

50 

10, 7 




-116 

-127 

10, 8 




-54 

3 

10, 9 




55 

70 

10, 10 




2 

3 

11, 0 





210 

n, l 





44 

n, 2 



1 


206 

n, 3 





0 

11, 4 





-67 

n, 5 





-30 

11, 6 





28 

11, 7 





2 

11, 8 





47 

11 , 9 





-11 

11, 10 





-34 

11, 11 





1 

12, 0 





-61 

12, 1 





162 

12, 2 





-228 

12 , 3 





75 

12, 4 





81 

12 , 5 





-69 

12, 6 





-64 

12, 7 





-7 

12, 8 





31 

12, 9 





30 

12, 10 





14 

12, 11 





•2 

12, 12 





0 
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S ^ e ^ cal - Har «onic Analysis of the Geomagnetic Field for 1955, 

• S * S * R ‘ ' W ° r JJ (3 ^ ta \ Mean va i ues oflL a nd Y , for series terminating 
with indicated values in P m . * 

n 
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lteble 3(^)cont'd 
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Ttble 3(c) Spherical-Harmonic Analysis of the Geomagnetic Field for 1955, 
U.S.S.R. World Charts . Mean values of Z , for series terminating 
with indicated values in P®. 


n, m 

Gaussian (Schmidt) Coefficients, g m (Units of 10”^ cgs) 

n 

a 

A 

A 

P 10 

10 

p 12 

12 

0, 0 

1196 

1772 

1733 

1751 

1740 

1, 0 

-301699 

-304698 

-304281 

-304319 

-304288 

1, 1 

-22973 

-22099 

-22126 

-22117 

-22117 

2, 0 

-15459 

-14476 

-14572 

-14530 

-14546 

2, 1 

29565 

29863 

29822 

29817 

29817 

2, 2 

16357 

16431 

16434 

16433 

16434 

3, 0 

13345 

9 886 

10494 

10420 

10459 

3, 1 

-20485 

-18114 

-18180 

-18138 

-18138 

3, 2 

11724 

12225 

12190 

12213 

12211 

3, 3 

8849 

8747 


8755 

8754 

4, 0 

9651 

10841 

10693 

10756 

10739 

4, 1 

7185 

7699 

7621 

7564 

7568 

1*, 2 

5476 

5784 

5793 


5788 

4, 3 

-3204 

-3876 


-3898 

-3894 

4, l* 

3283 

3226 

3225 

3220 

3222 

5, 0 


-3761 

-2969 

-3074 

-3019 

5, 1 


3472 

3372 

3463 

3466 

5, 2 


1355 

1271 

1320 

1295 

5, 3 


-546 

-409 

-462 

-467 

5, 4 


-1204 

-1293 

-1280 

-1278 

5, 5 


-1016 

-1058 

-1060 

-1058 

6, 0 


1380 

1194 

1281 

1259 

6, 1 


685 

562 

438 

46o 

6, 2 


616 

633 

622 

617 

6, 3 


-2427 

-2466 

-2489 

-2465 

6, 4 


-376 

-382 

-410 

-402 

6, 5 


358 

378 

386 

386 

6, 6 


-723 

-725 

-725 

-725 

7, 0 



905 

774 

846 

7, 1 



-135 

8 

13 

7, 2 



-137 

-51 

-108 

7, 3 



364 

250 

239 

7, 4 




-343 

-337 

7, 5 



-325 

-346 

-332 

7, 6 



195 

190 

192 

7, 7 



298 

301 

302 

8, 0 



-213 

-106 

-136 

8# 1 



-166 

-362 

-320 

8, 2 



26 

9 

-4 

8, 3 



-79 

-123 

-70 

8, V 



-22 

-92 

-74 


T 


,k 
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Ifcbl* 3(c)cont'd 


n, m 

Gaussian (Schmidt) Coefficients, (Units of 10 cgs) 

< 

4 

*8 

P 10 

10 

P 12 

12 

8, 5 



no 

166 

164 

8, 6 



-l6 

-25 

-14 

8, 7 



10 

-1 

-2 

8, 8 



182 

176 

176 

9 , 0 




-148 

-59 

9, 1 




194 

202 

9 , 2 




129 

30 

9, 3 




-207 

-226 

9 , 4 




142 

154 

9 , 5 




-91 

-47 

9, 6 




-30 

-5 

9 , 7 




36 

48 

9, 8 




-149 

-146 

9 , 9 




49 

43 

10, 0 




122 

86 

10, 1 




-259 

-196 

10, 2 




-26 

-48 

10, 3 




-73 

22 

10, 4 


i 


-149 

-n4 

10, 5 




168 

163 

10, 6 




-45 

-5 

10, 7 




-76 

-87 

10, 8 




-6 7 

-58 

10, 9 




38 

36 

10, 10 




19 

22 

11, 0 





101 

11, 1 





10 

11, 2 





-145 

n, 3 





-31 

11, 4 





24 

n, 5 





no 

11, 6 





89 

11, 7 





72 

11, 8 





23 

11, 9 





-68 

11, 10 





76 

li, n 





-33 

12, 0 





-43 

12, 1 





79 

12, 2 





-33 

12, 3 





147 

12 , 4 





64 

12, 5 





-11 

12, 6 





118 

12, 7 





-47 

12, 8 





64 

12, 9 





-14 

12, 10 





50 

12, 11 





14 

12, 12 





1 


’u’sVru 3 ^ e £ cal - 1 , «'”°nl= An.lyei» or the Geomghetlc Held for 1955 
U.S.S^R World Cherts Meen values of z , ^ eerie. tenstotiS 
with indicated values in P m 
n’ 

Gaussian (Schmidt) Coefficients, h“ (Units of 10 " 6 egs) ~ 


0 

0 

57984 

0 

-I8287 

3315 

0 

-5236 

2179 

-1333 

0 

1525 

-3749 

-846 

-1608 


0 

0 

58093 

0 

-18511 

3473 

0 

-4940 

2543 

-1337 

0 

1137 

-3086 

-857 

-1630 

0 

434 

982 

-20 

-1511 

1104 

0 

-516 

1327 

-38 

-146 

-81 

-432 


0 

0 

58174 

0 

-18590 

3445 

0 

-4740 

2553 

-1331 

0 

986 

-3176 

-914 

-1631 

0 

735 

1006 

21 

-1522 

1113 

0 

-756 

1166 

-224 

-154 

-30 

-423 

0 

408 

38 

109 

-48 

73 

-149 

-554 

0 

-323 

-250 

-377 

-PQ 


0 

0 

58170 

0 

-18594 

3454 

0 

-4759 

2524 

-1321 

0 

943 
-3139 
-895 

-1629 

0 

694 

944 
60 

-1520 

1114 

0 

-848 

1251 

-181 

-145 

-30 

-422 

0 

344 

-69 

194 

-38 

91 

-116 

-546 

0 

-468 

-107 

-293 


0 

58164 

0 

-18593 

3455 

0 

-4752 

2522 

-1329 

0 

938 

-3139 

-897 

-1632 

0 

740 

926 

28 

-1511 

1112 

0 

869 
1244 
-189 
-156 
•2 6 
-422 
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Table 4 
-2 

Weighted rms Error E in Units of 10 cgs in the Fit of Tabular Values 
at 10° Intervals of Co latitude and Longitude; U.S. and U.S.S.R. Charts 
for 1955 for Spherical-Harmonic Series Terminating in Degree 4, 6, 8, 
10, and 12 for (a) U.S. Charts and (b) U.S.S.R. Charts 



Component of Field 


Degree 

X 

Y 

z 

Mean 

(a) 

4 

0.223 

0.222 

0.324 

0.256 

(b) 

0.241 

0.295 

0.318 

0.285 

(a) 

6 

0.108 

0.116 

0.107 

0.110 

(b) 

0.151 

0.226 

0.186 

0.188 

(a) 

8 

0.098 

0.092 

0.071 

0.087 

(b) 

0.134 

0.220 

0.148 

0.167 

(a) 

10 

0.091 

0.097 

0.058 

0.082 

(b) 

0-132 

0.223 

0.129 

0.161 

(a) 

12 

0.087 

0.093 

0.054 

0.078 

(b> i 

0.131 

0.226 

0.101 

0.153 
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Table 5 

Values of Eqs . (16) and (18) for A0 = .04 and x q = 0 


(16) 

(18) 

5.02655 

4.44741 

.62832 

.60657 

.18617 

.18300 

.07854 

.07772 

.04021 

.03992 

.02327 

.02315 

.01465 

.01459 

.00982 

.00978 

.00690 

.00688 

.00503 

.00501 
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Table 6 


dS 

Values of J* for Various Heights Above the Sphere 
S 


z - 1 
o 

Value 

Center 

Error 

% 

Off Center 

Error 

% 

10. 

.00952 

.00948 

.4 

.00948 

.4 

l. 

2.09440 

2.10289 

- .4 

2.10312 

- .4 

.1 

54.39988 

54.02121 

.7 

53.99732 

.7 

.01 

619.00356 

626.65094 

-1.2 

705.91803 

-14. 

» 



Position Coordinates Estimated by Three Methods 
Field Line Starts at 40 N, 90 E; h * 0.1 Earth Radius 
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Table 8 

Components of Geomagnetic Field Estimated by Three Methods 
Field Line Starts at 40 N, 90°E; h - 0.1 Earth Radius 
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